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We present further Spitzer Space Telescope observations of the recurrent nova 
RS Ophiuchi, obtained over the period 208-430 days after the 2006 eruption. 
The later Spitzer IRS data show that the line emission and free-free continuum 
emission reported earlier is declining, revealing incontrovertible evidence for the 
presence of silicate emission features at 9.7 and 18 jttm. We conclude that the 
silicate dust survives the hard radiation impulse and shock blast wave from the 
eruption. The existence of the extant dust may have significant implications for 
understanding the propagation of shocks through the red giant wind and likely 
wind geometry. 

Subject headings: stars: individual (RS Oph) — novae, cataclysmic variables - 
binaries: symbiotic — binaries: close — infrared: stars 



1. Introduction 



RS Oph is a recurrent nova (RN) that erupted in 1898, 1933, 1958, 1967, 1985 flWallerstein 



20081 ). It consists of a semi-detached binary (orbital period 455.7 days) c omprising a roche 
lobe-filling red giant (M2III) and a massive ( > 1.2 M ) white dwarf fWD: lFekel et al.lboooh 
The eruption follows a thermonuclear runaway on the surface of the WD (jstarrfield et al. 
19881 ). In the case of the RS Oph cla ss of RN, however, the ejected material r uns into, and 



shocks, a dense red giant (RG) wind (IBode Kahnlll985l ; lO'Brien et al.lll992l ). 



The 1985 eruption was, for the first ti me, the subje ct of a mult i- wavelength observational 
campaign, from the radio to the X -ray (jBodd 119871 ). Its most recent eruption, on 2006 
February 12.83 (INarumi et al.ll2006l we take this to define the origin of time post-outburst), 



was the su 
(Das et al. 



j ject of an even more inten sive observational campaign. Infrared (IR) observations 
2007aJJb) showed evi dence for the shock, seen also at radio 



2006 



(lEvres et a 



Ness et al. 



2007 



2007 



Evans et al 



O'Brien et a 



Osborne et al 



2006) and X-ray (Bode et all 1200ft ISokoloski et all 12006 



20071 ) wavelengths, as the ejecta interacted with the RG 



wind. 



Observations of RS Oph obtained with the Spitzer Space Telescope ( Spitzer AWemei et al. 



20031 : iGehrz et al.l 120071 ) during t he period from 2006 April 16 through 26 UT (~ 67 days 



after outburst) were described by lEvans et al.l (l2007bl ). We present here further IR observa- 
tions conducted as part of a long duration sy noptic campaign, ob tained later in the outburst 
with the Spitzer Infrared Spectrometer (IRS: iHouck et al.l 120041 ) . 
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2. Observations 

RS Oph was observed with the Spitzer IRS as part of the Director's Discretionary Time 
program, Program Identification (PID) 270 on 2006 September 9.83 UT (day 208), 2006 
October 10.4 UT (day 246) and 2007 April 19.5 UT (day 367.2). Observations were performed 
using all IRS modules and the blue peak-up on RS Oph. The spectroscopy consisted of 5 
cycles of 14 second ramps in short-low mode, 5 cycles of 30 second ramps in both short-high 
and long-low modes, and 5 cycles of 60 second ramps in long-high mode. IRS basic calibrated 
data products (BCDs) were processed with versions 14.4.0 (2006 observations) and 16.1.0 
(2007 observation) of the pipeline. The extracted spectra w ere not defringed. Full details of 



the data reduction process are given in Evans et al.l (j2007bl ) and are not repeated here. 



Fig. [T] shows, with the same flux d e nsity s cale and over the wavelength range 6-14 /mi, 



the Spitzer spectrum from lEvans et al.l (l2007bl ). when free- free radiation and fine structure 



lines dominated the emission, along with the three new spectra reported in this paper. In the 
period covered by the observations reported here there seems to have been little change in 
the level of the emission: the differences of ~ 10% are comparable with the uncertainty in the 
flux calibration. We conclude that, in 2006 April, emission from the hot gas dominated the 
dust emission but as the emission from the gas subsided, emission from dust was revealed. 
The low resolution Spitzer spectrum of 2006 Sept. 09.83 UT is shown in Fig. 2, in which the 
silicate features at 9.7 /zm and 18 /zm are clearly visible. 

We confine our discussion in this paper to the silicate features, and defer a discussion 
of the evolving line and (gas) continuum emission to a later paper. 



Discussion 



The forma tion of silicate dust requires an environment in which O > C by number 
(IWhittetl 120031 ). Therefore on the assumption that the silicate dust we see in the Spitzer 
spectrum arises in the RG wind (see below), we conclude that the RG in the RS Oph 
system cannot be carbo n-rich; this is consistent with the presence of TiO bands in the 
optical (IBlair et al.lll983l ). and the presence of d eep O I absorption combined with no carbon 
absorption in the X-ray range (INess et al.ll2007l ). 



RS Oph is listed as a 12 /im source in the IRAS Point Source Catalog, with a flux density 
(not colour corrected) 0.42 Jy and upper limits in the other IRAS bands. To obtain a better 
estimate of the mid- IR flux densities, we retrieved the original scans from the IRAS database 
(lAssendorp et al.lll995l ) and used the GIPSY software to measure flux densities from traces 
through the position of the star. In this way the 12 /im value is refined to 0.37 ± 0.02 Jy, 



-4- 



and a marginal detection, 0.09 ±0.03 Jy, is recovered at 25 /im; RS Oph was not detected at 
60 or 100 /on, to uninter estingly large upper li mits. These data, alon g with the 10 /im flux 
from Geisel et al. ( 1970 ). the 2MASS survejll] ( Skrutskie et al. 2006 . data obtained on 22 
April 1999) and the Spitzer spectrum from 2006 September, are shown in Fig. [31 



Given that the IRAS, 2MASS and I Geisel et al.l (jl970j) data were obtained between 
eruptions we are confident that they represent the quiescent state of RS Oph. Moreover the 
fact that the Spitzer data are superficially consistent with the inter-outburst IR observations 
(see Fig. [3]), and that the dust was present at day 430, suggests th at the dust seen by Spitzer 
is likely present between outbursts (note that iHodge et al.l (120041 ) found no evidence of dust 
in ISO spectra of RS Oph, but their data only w ent as far as 9 /mi). We have made a 
preliminary fit using dusty ( llvezic fc Elitzurlll999l ) with a 3600 K blackbody (to represent 
the RG during q uiescence) as input radiation field. Using silicate optical constants from 
Ossenkopfi (119921 ). we find optical depth ry — 0.1 (± ~ 10%), dust temperature at the inner 
edge of the dust shell To ~ 600 K (± ~ 100 K) and dust shell inner radius of 1.5 x 10 14 cm. 
We derive a luminosity L* ~ 440 L^, mass-lo ss rate M ~ 2.3 x 1O~ 8 M yr _1 (typical for 
these stellar parameters; Ivan Loon et al.ll2005l ) for a dust to gas mass ratio of 0.005, and a 
terminal wind speed = 8.1kms _1 . 

While we caution against taking these derived parameters too literally, the 'outer' edge 
of the wind (defined by v^At, where At is the time between the 1985 and 2006 eruptions) 
is ~ 5.4 x 10 14 cm. The ejecta from the 2006 eruption, travelling at constant velocity, 
would have overtaken the outer edge of the wind even a t the earliest of our observations 
(t = 208 days). Ho wever the shock, which is decelerating (IBode et al.ll2006l ; iDas et al.ll2006 



Evans et al.ll2007af) may not, as yet, have r eached the wind edge. Taking Tv/V 9 . 7 ~ 17.8 
(cf. Mathis 19981 ) and silicate opacity from Ossenkopf ( 1992 ). we deduce a dust mass of 
4.1 x 10~ 9 M Q . Also, the inner radius of the dust shell is rather greater than the binary 
separation (~ 2.2 x 10 13 cm, assuming that the mass of the RS Oph system is ~ 2M ). 

We have determined the full-width at half-maximum (FWHM) of the 9.7/im feature by 
drawing a linear continuum between 8.5 /im and 12.5 /im (see Fig. HI in which the fluxes at 
10 /mi have been normalized to the 2006 September value). We find that, on 2006 September 
9 and October 17, the FWHM was 2.17/xm and 2.14 /im respectively, with no evidence for 
change; indeed it is difficult to distinguish the silicate features for these two epochs. By 2007 
April 19, however, the FWHM was significantly smaller, 1.99 /im, perhaps because the dust 
had experienced some processing (resulting in a change in the dust optical properties) since 



1 The 2MASS survey is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by NASA and the NSF. 
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2006 September/October. The profile of the 9.7 /im fea ture in RS O p h bro adly resembles 
those of the "Type C" silicate features of AGB stars in ISpeck et all s (120001 ) classification. 
However the silicate feature in RS Oph may dffer in detail from those of AGB stars for 
several reasons, such as lack of hot dust in RS Oph, or to grain size, shape, and composition 
differences. Differences may also arise because the dust in RS Oph is persistently exposed 
to ultra-violet (UV) radiation from the WD, as well as suffering frequ ent exposure to ho t 
ionized gas (see below); the latter is known to anneal amorphous silicate (ICarrez et al.ll2002l ). 



Dust formation in the winds of classical novae during eruption is well-documented 
(jEvans &: Rawlingd 120081 ; iGehra 120081 ). Indeed our observations were obtained during an 
inflection in the visible light curve, where the \^-band magnitude abruptly decreased by 
~ 0.6 mag, recovering shortly thereafter. In classical novae such inflections, accompanied by 
changes in the IR spectral energy distribution, are associated with dust condensation in the 
ejecta. 

However, in RS Oph we are confident that dust did not form in the ejecta , for two 
reasons. First, the temperature of the shocked gas (~ 10 6 K; lEvans et alJ l2007afl ) is not 
conducive to grain formation, which requires dense, relatively cool (~ 2000 K) gas. Second, 
the strength of the 18 /mi feature relative to that of the 9.7/zm feature indicates that the 
dust has been processed. In freshly-condensed silic ate dust, the 18^m feature is expected 
to be weak by comparison with the 9.7 /xm feature (INuth &: Hechtlll990l ). Alternatively, the 
18/im/9.7/im ratio for silicate features in these circumstellar environments could be solely 
dependent on dust temperature. We have already noted the apparent paucity of hot dust 
in RS Oph, so the 9.7/im feature will be less strong. In fact our dusty fit (Fig. [3]) over- 
estimates the 9.7 //m feature strength, even though the model satisfactorily fits the dust 
continuum as well as the 18 fim feature. We conclude that the silicate emission evident in 
the RS Oph Spitzer spectra is inconsistent with dust formation in the 2006 eruption. 

Previous mode l s of the evolution o f the remnant of RS Oph following eruption (e.g. 
Bode fc Kahnl Il985t lO'Brien et al.l Il992l ) have, for reasons of simplicity, assumed spherical 
symmetry. Accordingly the ejecta sweep into the RG wind, are shocked, and eventually reach 
the outer edge of the wind. In this scenario the entire RG wind is completely engulfed by the 
ejecta, and is replenished only after the eruption has subsided and the RG wind starts anew. 
As a result any dust condensing in the RG wind would also be engulfed by the shocked gas. 

However the presence of warm silicate dust so soon after the 2006 eruption, and its 
likely presence between eruptions, may imply that a substantial portion of the RG wind 
never 'sees' the eruption at all. Assuming that the RG fills its Roche lobe, its radius is ~ 0.7 
of the binary separation for a mass ratio (secondary/primary) ~ 0.4. It therefore presents a 
solid angle ~ 0.3 steradians to the WD (i.e. it covers ~ 2.5% of the WD's sky). 



- 6- 



When the RN eruption occurs on the surface of the WD it sweeps past the RG in ~ 1.5 
days, potentially arriving at the nominal outer edge of the wind in a few tens of days. This is 
<C than the orbital period, so ~ 2.5% of the wind is shielded by the RG and sees neither the 
eruption itself nor the expanding ejecta. However the situation is more complex than this 
in view of the radial- dependence of the wind density, and it is likely that the shocked gas 
will encroach into the shielded volume. Nevertheless our estimate of 2.5% for the fraction 
of shielded wind is likely of the right order: shielding of the wind by the RG is unlikely to 
account for the dust we see. 

On the other hand, the density of circu mstellar mat e rial in binaries is greatest in the 



binary plane (e.g. ISpruit fc Taamll200ll . and lGehrz et al.l (120011 ) for the case of the massive 
binary RY Set), so the RG wind will not be spherically symmetric. Outflow, both of the 
wind and the 2006 ejecta, will be inh ibited in the binary plane and in deed, anisotropies are 



suggested by H ST (IBode et a" 



interferometry (IMonnier et al 



20071) and radio (jO'Brien et al.l 120061 ) observations, and IR 



20061 ). We expect a significant amount of dust to survive in 



the binary plane; s uch a dense medi um woul d be consisent w i th the shock decele ration seen 
in the early Swift (IBode et al.ll2006l ) and IR ( IDas et al.ll2006l ; Evans et al.ll2007al ) data. 



The higher density material in the binary plane is effectively shielded and may not 
experience the eruption (which in any case may not have been spherically symmetric); this 
offers, in part, an explanation for the continued presence of silicate dust in the environment 
of RS Oph after eruption. 

Much of the dust in the RG wind away from the binary plane would not survive the 
surge of UV radiation emitted at outburst. The temperature of grains at the inner edge of 
the dust shell would rise by a factor [LEdd/-^*] 1 ^ 5 — 2.5 when heated by UV radiation from a 
source at the Eddington limit, L^d, for a Chandrasekhar-mass WD. Thus grains at 600 K at 
quiescence would be heated t o ~ 1500 K, in excess of the nominal evaporation temperature 
(~ 1300 K; ISpeck et al.ll2000l ) of silicate dust. Indeed, silicate grains would not survive the 
eruption out to distance r eva p ~ 2.2 x 10 14 cm from the WD in outburst. Furthermore, it 
is likely that a substantial amount of dust may have survived in denser regions swept-up in 
winds from previous outbursts. 

Any dust beyond the evaporation distance would survive the UV blast, so we con- 
sider w hether this dust wo uld survive the passage of the shock. As discussed in earlier 
papers fevans et al.l l2007al lbh. the gas temperature implied by the IR fine structure and 
coronal lines likely ranges from ~ 150 000 K to ~ 900 000 K, with an electron density 
~ 2.2 x 10 5 cm" 3 in the cooler re gion; the x-ray data indicate shock temperatures even 



higher than this (IBode et al.ll2006l ). Dust grains in this environment will be sputtered as 
they are engulfed by the shock. Assuming that the sputtering threshold for silicate grains is 
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~ 6 eV (ITielens et al.lll994l ). the time-scale for complete erosion of 0.1 fim silicate grains in 
the environment implied by the IR data is ~ years, where the yield Y (the number of 
grain atoms ejected per incident ion) is ~ 10~ 3 (ITielens et al.lll994j ). As this erosive time- 
scale is 3> At (the inter-outburst interval) any grains in the RG wind would survive passage 
of the shock, but might well be processed during the passage of the hot gas. 

Our conclusion is consistent with iV-band interferometric observations conducted on 



day 3.8, which may have independently seen the dust we report here ( iBarry et al.l 120071 ). 
The nearest constructive peak in the interferometer's acceptance pattern is at an angular 
distance from the central null fringe that depends on wavelength as X/B, where B is the 
length of the projected baseline (84 m). At the center of the iV-band, in which Si I 9.407/zm 
and dust are detected, the peak transmission is at angular separation ~ 12.5 mas from the 
center of the source brightness distribution, equivalent to ~ 3 x 10 14 cm at a distance of 
1.6 kpc. The initial shock velocity would need to have been ~ 9000kms _1 to reach this 



distance, travelli ng at constant velocity, far higher than observed (e.g. lEvans et al.ll2007al ; 
Bode et al.ll2007l ). Therefore this emission line was not excited by passage of the shock front 
but by the UV flash from the TNR event. This measurement is consistent with our assertion 
that dust detected by Spitzer was not created by the collision of the ejecta with the RG wind 
or by subsequent condensation, but was more likely formed during quiescence. 

Our continuing Spitzer observations of RS Oph will enable us to examine the effects of 
irradiation on the dust, and of the re- formation of dust as the circumbinary environment is 
re-established by RG wind infill. 



4. Concluding remarks 

We have presented further Spitzer IRS observations of RS Oph in the aftermath of its 
2006 eruption. Unlike our earlier observation, which showed emission by the hot shocked 
gas, more recent IR spectra reveal the presence of silicate dust. 

The dust we see could not have been formed in the 2006 eruption. It most probably 
survived the 2006 eruption, because a portion of the RG wind (including the dust) did not 
see the eruption at all, and because the passing shock failed to destroy any dust that survived 
the UV blast at the eruption. This conclusion may have major implications for the evolution 
of the shock, and of the long-term survival of the RG wind between eruptions. We predict 
that the dust we see around RS Oph will persist until the next eruption. 

This work is based on observations made with the Spitzer Space Telescope, which is 
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Fig. 1. — Spitzer IRS spectra of RS Oph over the period 2006 April (hot gas dominant) to 
late 2006 and later (silicate features prominent). 
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Fig. 2. — Spitzer low resolution IRS spectrum of RS Oph obtained on 2006 September 
9.8 UT. The silicate features at 9.7 and 18 /xm are clearly visible. 
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Fig. 3. — Spitzer I RS spectrum from 2006 September 9.8 UT, together with data from 
2MASS, IRAS, and lOeisel et al I Jl97oh . Triangles: 2MASS; squares with error bars: IRAS; 
large open square: iGeisel et al.l (119701 1 . The broken curve is a 3600 K blackbody, used as 
input into DUSTY; the dotted curve is the dusty fit. For details see §3. 
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Fig. 4. — Spitzer spectrum of RS Oph for 2006-2007; a linear continuum has been removed 
from the 8.5-12.5 /im region. Principal H recombination lines (n — m) are identified. Spectra 
for 2006 September and October are essentially identical; note that the 9.7 //m feature for 
2007 April (dotted line) is noticeably narrower. 



